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Abstract 

We study the scattering of two Skyrmions at low energy and large separation. 

O ■ We use the method proposed by Manton for truncating the degrees of freedom 

O I of the system from infinite to a manageable finite number. This corresponds to 

Q\ ', identifying the manifold consisting of the union of the low energy critical points 

of the potential along with the gradient flow curves joining these together and 

A ■ by positing that the dynamics is restricted here. The kinetic energy provides 

an induced metric on this manifold while restricting the full potential energy to 

the manifold defines a potential. The low energy dynamics is now constrained 

K> , to these finite number of degrees of freedom. For large separation of the two 

j^ ■ Skyrmions the manifold is parametrised by the variables of the product ansatz. 

We find the interaction between two Skyrmions coming from the induced metric, 

which was independently found by Schroers. We find that the static potential is 

actually negligible in comparison to this interaction. Thus to lowest order, at large 

separation, the dynamics reduces to geodesic motion on the manifold. We consider 

the scattering to first order in the interaction using the perturbative method of 

Lagrange and find that the dynamics in the no spin or charge exchange sector 

reduces to the Kepler problem. 






1. Introduction and Conclusion 



Scattering of solitons in a non-integrable, non-linear classical or quantum field 
theory remains an intractable and difficult problem. It, however, concerns one of 
the most interesting aspects of the nature of the corresponding physics. Numerical 
methods have given reasonable ideas on how the scattering proceeds but they are 
still unsatisfactory for uncovering the detailed dynamics governing the scattering. 

A method has been proposed by Manton^ for truncating the degrees of free- 
dom from the original infinite number to a relevant finite number of variables for 
the case of low energy scattering. One first considers theories of the Bogomolnyi 
type. In these theories there is typically a topological charge (an integer A^) which 
caracterizes the soliton sector giving in some sense the number of single solitons 
in the sector. In each topological sector there exist exact static solutions with 
a fixed number of arbitrary parameters, the moduli space, which describe these 
solutions. The number of such parameters is equal to the number of parameters 
for the one soliton solution multiplied by the topological charge (the number of 
solitons). There exist asymptotic solutions which are easily identifiable with A^ 
single solitons with large mutual separations, but for small separations the indi- 
vidual solitons are subject to strong deformations and they lose their identity. As 
these are static solutions at essentially arbitrary separations between the solitons, 
clearly there can be no static forces between them. Otherwise the solutions would 
not exist; the solitons would move towards or away from one another. Thus in 
each soliton sector the sub-manifold of solutions corresponds to an equi-potential 
surface, the minimal energy surface. The kinetic energy can be properly inter- 
preted as corresponding to a metric on the space of all configurations. Typically 



we have 



d^x -gij{(f)k) ^i^j- (0) 



In many instances Qij is just flat, but it can be non-trivial, for example, for non- 
linear sigma models such as the Skyrme model. In any case, in general, the induced 
metric on the sub-manifold of static solutions is non-trivial. If one considers a time 
evolution with initial conditions corresponding to a point on the sub-manifold, with 
initial velocity tangent to the sub-manifold and arbitrarily small, it is intuitively 
evident that the evolution will remain on the sub-manifold and correspond simply 
to geodesic motion. 

It is a difficult task to prove that such a truncation of the degrees of freedom is 
a sensible, mathematically rigorous, perturbative scheme. The non-linearity of the 
theory means that the degrees of freedom tangential to the equi-potential surface 
are coupled to all other degrees of freedom at higher order. We are assuming that 
these couplings are negligible. Gibbons and Manton^ applied this program with 
remarkable success to the case of magnetic monopoles in the BPS limit and it has 
also been applied to many other situations, including vortex scattering, Kaluza 
Klein monopoles, black holes, lumps in CP^ models, etc. (We refer the reader to 
the article of Samols"^ for detailed references). 

The generalization to the more common situation where the set of static so- 
lutions do not include solitons at essentially arbitrary separations can also be 
developed. Here the forces between the solitons do not exactly cancel, but these 
are assumed to be weak. The moduli space of the minimum energy critical point 
is of a smaller dimension than before and does not ever represent many isolated 
single solitons. Thus truncation of the dynamics to the sub-manifold of such con- 
figurations is inadequate to describe the scattering of the solitons. It is, however, 
evident that the moduli space of the asymptotic critical point of A^ infinitely sep- 



arated single solitons will in fact he N x k dimensional, where k is the dimension 
of the moduli space for one soliton. It is argued by Manton, that the low energy 
scattering will be restricted to or will lie close to the part of the configuration 
space corresponding to the union of the gradient fiow curves which start at the 
asymptotic critical point and reach other critical points. They will reach the mini- 
mal critical point or, in fact, other critical points, which may actually have energy 
more or less than the starting point. We find this idea also intuitively reasonable. 
As long as the gradients involved are not too steep very little radiation should 
be incurred. For the Skyrnie model, numerically it is known^ in the two soliton 
sector that the difference in energy between the "deuteron" bound state and the 
asymptotic, infinitely separated, configuration is rather small, thus the gradients 
involved should not be large. 

In this paper we identify the sub-manifold corresponding to the gradient flow 
curves for the Skyrme model, in the two Skyrmion sector for large separation 
between the Skyrmions, as the set of configurations parametrized by the prod- 
uct ansatz (undeformed). The force between two Skyrmions at large distance 
is governed by their far field behaviour. For the product ansatz this gives rise 
to an interaction which falls off like 1/d^, where d is the separation. Localized 
deformations of each Skyrmion will not affect the behaviour of the interaction. 
Modification of the far field behaviour of each Skyrmion will only weaken their 
mutual interactions if the self energy of the modified Skyrmions is to be kept fi- 
nite. Modifications which sharpen the far field fall off are permitted while those 
which weaken the fall off cause the self energy to diverge. We are not concerned 
here with intermediate range modifications and interactions. Therefore it is clear 
that the product ansatz with undeformed Skyrmions at large separation gives a 
good parametrization of the manifold of gradient fiow curves up to order 1/d^. 
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This fact is confirmed by numerical calculations where it is observed that the un- 
deformed product ansatz is surprisingly accurate even down to separations of the 
order of the size of a Skyrmion"^. 

We compute the induced metric on this manifold coming from the kinetic 
energy. We actually find that the metric dominates the interactions between the 
Skyrmions, the static potential is relatively negligible. To first order in the inverse 
separation the free metric is modified by a term which behaves like 1/d, which is 
much greater than the potential which behaves like l/rf"^. The characteristic scale 
is set by /tt, by e and the details of the Skyrmion configuration. It is reasonable to 
expect that it is given by the size of the Skyrmion. We remark that this scale has 
nothing to do with the scale given by a light massive pion. We have considered 
a massless pion, adding a small mass does not drastically modify the Skyrmion 
size. At large distance of course the pion mass cuts off all interactions with the 
usual Yukawa exponential. Thus our results are valid for the separation in the 
range that is large compared with the Skyrmion size but small compared with the 
Compton wavelength of the pion. 

The scattering, to first order, is again described by geodesic motion on the 
manifold parametrized by the product ansatz variables. Even though this is a 
great simplification from an infinite number of degrees of freedom to just 12 col- 
lective coordinates, the problem remains intractable. We calculate 2 conserved 
"angular" momenta coming from invariance under right iso-rotation coupled with 
spatial rotation and single left iso-rotation. Then we use the Lagrange pertur- 
bative method to obtain approximate equations of motion which correspond to a 
systematic expansion in the inverse separation. It is actually quite reasonable to 
make this further approximation since we have already dropped terms of second 
order in the induced metric. The equations are still rather complicated and we 



expect a rich structure in the scattering cross section. We show that the system 
reduces simply to the Kepler problem for Skyrmions synchronously rotating in a 
direction orthogonal to the scattering plane. 



2. The Skyrme model 



The Skyrme model is described by the Lagrangean density, 



^ triU^d^UU^d^U) + ^^ 



'Csk = -^ tr{U^d^UU^d^U) + ,^^ tr{[U^d^U, U^d^f) (1) 



where U{x) is a unitary matrix valued field. We take 

U{x) G SU{2). 

The Skyrme Lagrangean corresponds to the first two terms of a systematic expan- 
sion in derivatives of the effective Lagrangean describing low energy interaction of 
pions. It should be derivable from QCD hence f-j^ and e are in principle calculable 
parameters. These calculations are actually unfeasable and we take f^^ and e from 
phenomenological fits. What is even more surprising is that the Skyrme model 
includes the baryons as well. These arise as topological soliton solutions of the 
equations of motion. The original proposal of this by Skyrme^ in the 60's was put 
on solid footing by Witten^ in the 80's. For further references and details see the 
review article by Wambach and Walhout^. 

The topological solitons, called Skyrmions, correspond to non-trivial map- 
pings of M^ plus the point at infinity into SU{2): 

U{x):]R^ + oo^SU{2) = S^. (2) 
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But 

R^ + cx) = ^^ (3) 

thus the homotopy classes of mappings 

U{x) :S^^S^ (4) 

which define 

UsiS"") = TL (5) 

characterize the space of configurations. 

The topological charge of each sector is given by 

TV = -^ [ d^xe'^'' tr{U^d,UU^djUU^dkU) (6) 

247r^ J 

which is integer and is identified with the baryon number. Thus for the scattering 
of two Skyrmions, we are looking at the sector of baryon number equal to 2. In this 
sector the minimum energy configuration should correspond to the bound state 
of two Skyrmions, which must represent the deuteron. The asymptotic critical 
point corresponds to two infinitely separated Skyrmions. There exist, known, non- 
minimal critical points, corresponding to a spherically symmetric configuration, 
the di-baryon solution^. The energy of this configuration is about three times the 
energy of a single Skyrmion. There are also, possibly, other non-minimal critical 
points with energy less than two infinitely separated Skyrmions^. The scattering 
of two Skyrmions will take place on the union of the paths of steepest descent 
which connect the various critical points. 



3. Lagrangean of the Skyrmion-Skyrmion system 



With the definitions 



^1{U) = -- tr[r-d^UU^] 






(7) 

(8) 
(9) 



the Skyrme Lagrangean density becomes 



J TT 



Csk = '-fC-{u)^-c{uY 



1 

4? 



c{u)^ ■ c{uY c{u), ■ c{uY - c{u)^ ■ c{ur c{u), ■ c{uy 



(10) 



where the ■ impUes a contraction in isospace. We can separate this into a kinetic 
energy T which is the part quadratic in time derivatives and a potential energy V 
without any time derivatives: 



T 



V 



f^ 1 

-£o-/:o-^ 



f 1 



£0 • Ci Co • Ci — C{] • Co Ci • Ci 



^i ' ^i ^j ' ^j ^i ' ^j ^j ' ^i 



(11) 



(12) 



We consider the scattering only for large separation, thus we do not have to 
know the structure of the complicated region where the two Skyrmions interact 
strongly and consequently are much deformed. In the region of large separation 
the product ansatz corresponds to 



U{x) =Ui{x-Ri) U2{x-R2) 

= AU{x - Ri)A^ BU{x - R2)B^ 



(13) 



where U{x — Ri) and U{x — R2) correspond to the field of a single Skyrmion 
solution centered at Ri and R2 respectively. This particular ansatz gives a 12 
dimensional configuration space. It can be seen that 

£^(t/it/2) = C^iU^) + D„,(t/i) £^(t/2) (14) 

so T and V each separate into three parts: one depending only on JJi, another 
depending only on U2 and a third function of both Ui and U2 which describes 
the interaction between Skyrniions. This last term has already been investigated 
for the potential part V, and gives a contribution of the order 1/d^ where d = 
ll-Ri — -R2II (For details see Ref. 9 and references therein). We shall concentrate 
here on the interaction part of the kinetic energy T given by: 

1 



2e2 



{Cl ■ Cl Cl Cj + 1 ^ 2) +2{Cl- Cl + 1 ^ 2)Cl ■ D ■ Cj 

+2 {Cj-C^ + 1^2)Cl-D.Cl + 4Cl-D.ClCl-D.CJ 
-{Cl.D.Cif-{C]-D-Clf-2Cl-C]Cl-Cl 



(15) 



2Cl-D.ClC]-D.Cl 



2 {Cl ■C] + 1^2){Cl-D.C^ + Cl-D- Cl) 



where £^ = £^(t/i), Cf^ = £^(t/2) and D = DabiJJi). Since the time dependance 
of Ui is contained in A and i?i, we find that 

C-liUi) = (5'''-Dab{AUix-Ri)A^)] C^oi^) - Dab{A) A\ C\{U{x-Ri)) (16) 

and correspondingly for U2, B and R2- 

It can be seen that {5''^ -Dah{AU{x-Ri)A^)) is of order l/r^ while C'i{U{x-Ri)) 

is of order 1/r^ at large distance due to the behaviour of F which falls off like 1/r^. 
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This actually implies below that the term containing £^ above is not important 
to the term of leading order in the kinetic energy. 

We now outline the computation of the kinetic energy Tint obtained after 
integrating Tint over all space. The procedure is similar to the one used in Ref. 9 
to calculate the interaction from the potential. We divide space into three regions 
I, II an III. II is a spherical region centered on the first Skyrmion, of a radius 
much smaller than d but large enough so that outside region II the asymptotic 
behaviour of F is valid. Region III is the same for the second Skyrmion. Region I 
is of course the remaining space and the field coming from both Skyrmions behave 
asymptotically here. We find the leading contribution to Tint behaves as 1/d and 
comes from the first term in Tint evaluated in region I. The contributions coming 
from the Skyrme term are of higher order due to the behaviour of £^. The leading 
contributions from region II and III are of order 1/d'^ and hence negligible. The 
interaction over region II is actually computed by extending the integrand to the 
whole of space. This is justified since it only modifies the contribution at higher 
order in 1/d. 

This simplifies the evaluation of the integral and we find 



:,nt = / d^xf^C^oiUl)Dab{Ui)C^oiU2)+0{l/d^) 

A . . ^''^ 

rac^jhd ^ar^^ 'Jldr^^ /^^J _ ^t^j\ JJJA^B) + 0(l/d^) 

d 



where A = 27iK^f^, F{r) ~ n/r^ at large r, 7^"(A) = 7^g(A) and d = d/d. We 
have used the relation 71°'{A) = Dab{A) C^(A). The free part of the kinetic energy 
is well known^ and we can finally write the Lagrangean L for the N = 2 sector of 
the Skyrme model to leading order of the dynamics of the variables of the product 
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ansatz: 



L = -2M + ^Md^ + 2A(£"(A) £"(A) + £"(S) £"(S)) 

(Jj 



(18) 



where 



'o 



x^k 



9r 



+ 



1 sin^F 

2^ r2 



sin^F 



r 



zap 

2 \ ^^ 



2n 



(19) 



is the mass of a Skyrmion and 

A = (c/tt)^ / r'^dr sin^ F 



1 + 



(e/.: 



1/2 



F'" + 



sin^F 



(20) 



is its inertia momentum. We have replaced Ri — i?2 by (i placing us in the center of 
mass reference frame and reducing the number of degrees of freedom of the system 
to 9. Equation (18) is the main result of our calculation. 

The metric can easily be obtained from this expression by choosing local 
coordinates on the product ansatz manifold and extracting the quadratic form 
relating their time derivatives. With the potential part absent from the Lagrangean 
of the system (to first approximation) the solution of the problem now resides in 
finding the geodesies of the metric on the product ansatz manifold. 

The result of equation (18) was also independently obtained by Schroers^'^. 
He found a leading contribution which behaves as 1/d and even calculated sub- 
leading spin-orbit coupling terms. The only other comparable calculation to our 
knowledge has been done by Walhout and Wambach^ for the case of massive 
pions. The limit as myr ^ of their expression, however, does not leave a term 
which behaves as 1/d and hence does not reproduce our result. We believe that 
this contribution should come also from their evaluation of the integral giving the 
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induced kinetic energy in the far field region (region I) and then recovering our 
result as m-j^ -^ 0. Such a contribution would also be proportional to (e""^'"'^)^, 
what they call "two pion exchange" . We believe that they have not computed the 
contribution from this region. 

We also add that our term is leading order in an expansion in inverse sepa- 
ration with respect to a scale which has nothing to do with the pion mass. This 
means that there are two length scales for nucleon-nucleon interaction predicted by 
the Skyrme model. It would be interesting to see if this can be phenomenologically 
or experimentally justified. 



4. Approximate Euler-Lagrange equations 



Unfortunately the two body Lagrangean has a very complicated structure. 
Simply finding the expression of the Euler-Lagrange equations is a long and tedious 
undertaking (let alone finding solutions to these equations) . We do not record the 
equations here. The Lagrangean possesses numerous symetries and corresponding 
conserved quantities. There are two interesting symetries apart from translational 
invariance giving conserved charges. 

From left iso-rotation 



A^CA 

(21) 
B^CB 
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where C is a constant SU{2) matrix we find the conserved "total isospin": 

Jl = 4A[C''{A) + C''{B)] 
A 



^dl 



''^'D^,{A){d'^ - d'd^)e^''''D,,{B)n\B) (22) 



The Lagrangean is also invariant under right iso-spin rotation coupled with a 

spatial rotation 

A^AC 

B ^BC (23) 

where C is a constant SU{2) matrix and Dab{C) is its representative for spin 1. 
This gives the further conserved "total angular momentum" : 

J^ = M{dx d)'" + 4A[n''{A)+n''{B)] 

A ^^^ak^JbdJ^d^^^^g^J_^^^J•^^^^^J^tB) (24) 

^^^ac^J6fc^c(^)(^^i _ d' d^ ) D ah{A^ B) 



+7 



These conserved quantities do not help us greatly in solving the equations for the 
geodesies. In fact it should be noted that the peculiar form of the Lagrangean 
(not unlike one describing a pair of coupled rigid bodies floating in free space) 
does not even allow the usual separation of the rotational movement into a global 
and relative part. In order to go further we must resort to a supplementary 
approximation. 

We use the perturbation method of Lagrange^^ familiar in celestial mechan- 
ics. This approximation scheme neglects the changes induced on the free canon- 
ical momenta by the interaction term. To begin with we need the free Pois- 
son brackets for the degrees of freedom of the system. We use the variables 
(J^rfi£"(A),£"(S),7^"(A),7^"(S)). £" is proportional to the isospin of a 
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Skyrmion and 7^" to its spin. We can easily compute the Poisson brackets for 
the free theory (the hmit d goes to infinity of L) : 



1-j 



{d\U^}^d 






(25) 



{£"(A),7^^(A)} = 
{C^{A),DUA)} = -^e^''D,M) 

{n-{A),DUA)} = ^e^^'D.^iA) 

where 11* is the conjugate momenta to (i*. Because of the symmetric nature of the 
free Hamiltonian, the same brackets are true if we replace A by S everywhere, and 
all the mixed brackets between A and B are zero. The free system in the center of 
mass reference frame possesses 5 vectorial conserved quantities: (i*, £"(A), C"'{B), 
TZ°'{A) and TZ°'{B). Let C* be one of those quantities. Its time derivative is given 

by 

at 
where H = Hfree + Hint. The approximation comes by using the free Poisson 
brackets to compute {C\ Hint}o- This is of course only correct to first order. 
Since C* is conserved in the free system, {C*, if/^eejo = ^^^ 

With the free Poisson brackets we obtain without difficulty the following coupled 
system of equations: 
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d'' + 



2A 



ikji I sik 



S-Jd" + 5'''d'+5'^d' - Wd'd' 



e""'e'^'^n%A)n'^{B)Dab{A^B) = 






2Md 



t'''^t^^'^n'^{B){5'^ -d'd^) 



dt 



iV'iB) 



2Md 



e>'^fnf{A)Dab{A^B) + e''''fDfb{A^B)n^{A) 



^^^c^jbd^c^A){5'^ -d'd^) 



e'"^fnf{B)Dah{A^B) + t^^fDaf{A^B)n'^{B) 



(26) 



Our approximation is reliable as long as the separation d between the particles 
is large enough for the conjugate momenta to stay close to their free values. As 
we have already worked with the undeformed product ansatz approximation, and 
neglected the potential, which are both valid for large d, we feel confident that we 
have not lost any meaningful information by making this further approximation. 
If d is kept large we should then find geodesies similar (qualitatively at least) to 
those given by the exact equations of motion. The system of equations (28) is 
still quite complicated but some specialized solutions are easy to obtain and give 
simple trajectories as we shall now see. 



5. Skyrmion-Skyrmion scattering 



The simplest geodesic is obtained by taking B = A and replacing it in the 
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equations (28). We then get the foUowing system of equations: 

•* 2A 
d" + 



M(P 



d''(n''{A)n''{A) + 3(d"7^"(A)) j - 2n^{A)d"-Tr{A) 

j^C\A) = (27) 

^n\A) = — ^n''(A)d''e''^'d^n^(A) 

dt ^ ' 2Md ^ ^ ^ ^ 

We see that if A is chosen at t = such that d is paraUel or perpendicular to 
TZ^ (A) then the last two equations become C^ = and TZ'^ = and are satisfied if 
C'' (A) and TZ^ (A) are constants. This means that A corresponds to an iso-rotation 
about a fixed axis with constant angular velocity. 

Let us begin by choosing 71^{A) parallel to d. This does not actually lead 
to a scattering because d does not have the liberty to change direction with time 
(it has to stay parallel to 7^^ (A) which is a constant). Substituting the condition 
TZ^ = a d^ = constant in the equation for d we find an equation for a particle 
constrained on a line with an attractive Coulomb potential. The trajectory will 
necessarily lead to d small and hence to regions where we can no longer trust our 
equations or their predictions. 

The case where 7Z^{A) is perpendicular to d is more interesting. Again the 
last two equations are satisfied with A corresponding to a steady iso-rotation, 
however d now has the freedom to move in a plane according to the equation: 

d'' + ^^n''{A)n''{A)d^ = 0. (28) 

This is the equation of a Kepler system: the two Skyrmions scatter in a plane 
while keeping their isospins and spins constant and perpendicular to the plane of 
the orbit. 

These geodesies are the simplest to compute algrebraically but are almost 
all that we can compute by hand. Indeed if we take TZ'^ at t = not exactly 
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parallel or perpendicular to cl", complicated precessions appear. It seems that the 
typical motions of the Skyrmions described by equations (28) are very complex 
and allow for large nutations, precessions and spin flips. Quantum mechanically 
these would correspond to neutral and charged pion exchange which we expect 
from a dynamical model representing nucleon-nucleon interactions. 
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